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Abstract 

We construct a family of measures for random fields based on the 
iterated subdivision of simple geometric shapes (triangles, squares, tetra- 
hedrons) into a finite number of similar shapes. The intent is to con- 
struct continuum limits of scale invariant scalar field theories, by imitat- 
ing Wiener's construction of the measure on the space of functions of one 
variable. These are Gaussian measures, except for one example of a non- 
Gaussian fixed point for the Ising model on a fractal. In the continuum 
limits what we construct have correlation functions that vary as a power 
of distance. In most cases this is a positive power (as for the Wiener mea- 
sure) but we also find a few examples with negative exponent. In all cases 
the exponent is an irrational number, which depends on the particular 
subdivision scheme used. This suggests that the continuum limits corre- 
sponds to quantum field theories (random fields) on spaces of fractional 
dimension. 



1 Introduction 

Wiener constructed a Gaussian measure on functions of one variable by sub- 
dividing the interval and imposing the scale invariance of the measure [T]. By 
subdividing the plane into triangles, and space into tetrahedrons, we construct 
analogous measures for functions of several variables. These are field theoretic 
analogues of the hierarchical models [2] commonly studied in spin systems [3]. 
They can also be viewed as the quantum analogues of the Laplace equation on 
fractals studied by Strichartz|4J. The method we use is to seek fixed points 
of the "real space renormalization group" , which has been in use in statistical 
physics for many years [5]. 



* arnabkar@pas.rochester.edu 

trajeev@pas.rochester.edu 

' Also at the Department of Mathematics 



There are other approaches to constructing quantum field theories in frac- 
tional dimensions [S] , in part motivated by the success of dimensional regular- 
ization. But our construction appears to be more direct, being based on a finite 
element method. It does not appear to be possible to construct scale invariant 
Gaussian measures of quantum field theories in integer space dimension greater 
than one by our method of iterating finite integrals. We extend our ideas of scale 
invariance to other systems like the Ising model, to find a non trivial unstable 
fixed point for a fractal of dimension greater than two. The critical exponents 
are calculated in this case using ideas of statistical mechanics. 

In the continuum limit what we construct is a quantum field theory of frac- 
tional scaling dimension, with a correlation function that varies as a power of 
distance; this power is often positive, but we also find a few cases where it is 
negative. An example of a system with growing correlation function would be 
that of Brownian motion. The average distance two points undergoing Brownian 
motion increases with time as < [x(ti) — xfo)} 2 >~ \ti — t 2 \- 

Physicists denote the Wiener measure by 



Exactly what does this mean? There certainly is no Lebesgue measure which 
can be called T>4> in the space of functions because that space is infinite di- 
mensional. Also, the functions <p are, with probability one under the Wiener 
measure, not differentiable. cj)' is at best a distribution. So it does not have a 
sensible product: <fi' 2 is by itself meaningless. 

Recall that the derivative ^ in calculus is to be understood as a limit con- 

ax 

structed using epsilons and deltas of modern analysis, and not as some infinites- 
imal divided by another. Similarly, the path integral is to be understood as a 
limit of finite dimensional integrals. The action being the integral of a local La- 
grangian ( a function of the field and its derivative) is not to be taken literally. 
Indeed, a limit of a sequence of local measures (that depends only on nearest 
neighbors of the triangulation) may not tend even formally to a local action 
principle. 

In order to construct the scale invariant measures in higher dimensions, the 
field variable could acquire an anomalous dimension. As a consequence, the 
correlation functions have irrational numbers as their scaling exponent. In the 
simplest examples, G(xi,x 2 ) =< [4>{xi) — 4>(x 2 )] 2 > scales as | a^i — ;e 2 | - 611516 
for the subdivision of right angled triangled triangles and as |a;i — ;e 2 | - 329351 f or 
tetrahedrons. These exponents are not universal as other subdivision methods 
lead to other exponents We suggest that this exponent of the correlation function 
of the Gaussian be used as a measure of the dimension of a fractal, analogous 
to the Hausdorff dimension: it is more relevant for the behavior of randomness 
in quantum systems. 
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2 Subdividing the Real Line 



Suppose that cf>o and <f>\ are the boundary values of the field on the interval 
[xo, xi] of length a = x\ — xq. Integrating out the field everywhere on the inside 
of the interval induces some measure h(<j>i, Xi|</>o, Xo)d(f)od(j)\ on the boundary 
values. The simplest possibility is a Gaussian that is translation invariant in 
field values 

htfuxtlfaxo) = C oe -W*i-*>) a 

for some a a and C a which can depend on |xi — x |. If we subdivide the 
interval into two (for example, at its midpoint x i = \{x\ + x )) and integrate 
over the field at the new point, we should get back the original measure. 

f h(<f>i,xi\<f>oi,xoi)h(<t>oi,xoi\<f>o,%o)d<t>w = h{4>i, xi|«^o, %o) 
This leads to the conditions 



C| j exp|-^o;| [(<fc 



: 0io = C exp|--a (<^i 
C n = Cl 



<f>oY 



By completing the square, 



C\ 



I 



exp ■ 



->oi 



ho exp 



\2 a ^ 2 f = ^exp j--a a (>i - O ) j> 



Thus, 




Solving 



a a = 



a 
a 



C a 



27ra 



v /27r(xi - x ) 



-e 2 ( x i- a: o) 



(0i-0o) 2 
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By repeated subdivisions of the interval we would get the integral 



/ Y[d^e 2 ^ •*+!-* 

i 

In the limit of small size of the interval, the exponent can be interpreted as 

f / 

justifying the physicists' notation 

But we reiterate that, this is not to be taken literally: \\m Xi+1 ^ Xi f'**Zf' 

almost never exists and so 4> 12 is actually meaningless. The quantity a may be 
related to mass divided by h in some physical system or it may represent the 
inverse of the diffusion coefficient for a different system. 

A higher dimensional analogue of this process of subdivision would involve 
triangulating space-time. We would subdivide each triangle (or more generally 
simplex) to get smaller ones by introducing extra vertices. Integrating out the 
field at the extra vertices should give back the original measure. 



3 Subdividing Triangles 

Our idea is to approximate the functional integral of the quantum field theory 
by an integration of the scalar values at the vertices of a triangulation: 



/ V(h~ m A (<f>( X1 ), ■ ■ ■ 4{x n )) « / n n p - a Mao, • • • (j>(x n )) 

J J I f 

A triangulation is chosen with the points x a , a = 1, • • • n in the argument of 
A among the vertices. The measure of integration is assumed to be local; i.e., 
can be expressed as a product of factors Ff (4>o 7 4>i, ^2) each depending only on 
the field values at the vertices of a single face / of the triangulation. The basic 
idea is that if we subdivide a triangle by introducing new vertices, and integrate 
over the field at these new vertices, we must get back the contribution of the 
original triangle. 

What kind of triangle can be sub-divided into two equal triangles which are 
each similar to the original one? A moment's reflection shows that a right angled 
isosceles triangle has this property. Let us work this out in two dimensions to 
get a feel for the method before going to the much harder higher dimensional 
cases. However, it has been shown that for higher dimensions, one can subdivide 
a d simplex in n dimensions into n d similar simplexes. [7] 
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Figure 3.1: » 
We start with the triangle (Figure 3.1) with vertices xq = (0,0), x\ 



(a, 0), X2 — (a, a). Let 

be the factor corresponding to such a triangle, for some normalization con- 
stant C a and quadratic form Q a . We assume that it is invariant under transla- 
tions in the field variables for simplicity. It is likely that "massless" scalar fields 
may have this property. If the quadratic term was independent of cross terms 
between two vertices, it may explain scalar fields in the "infinite mass" limit. 
Also the interchange of the vertices 1, 2 must be a symmetry. Thus we have 

Qa(0O, 01, 0a) = "a ((00 ~ 0l) 2 + (02 ~ 0l) 2 ) + Pa{<t>0 ~ 2 ) 2 

for some constants a a , f3 a . This can also be expressed in matrix form as 

+ Pa -Ola ~Pa 
~Pa -OL a OL a + Pa 

If we choose the new internal vertices to be midpoints of the sides of the 
original triangle xqX\Xi 

1 . 

Xij = -(Xi + x 3 ), i,j =0,1,2 

it divides the original triangle into the four equal triangles that are similar to 
the original. The new triangles (a^oi^cs, ^02^01^1, ^1^12^02, £02^1222) have 
orthogonal edges which are half the length of the original triangle. 

The reproducibility condition obtained after integrating out the field vari- 
ables at the new vertices is 

2) _ / e -B, a {4>0,4>X,4>2, 001.002.012), 



5 



where, (up to a normalization constant) 



Ra(<h> 01, 02, 002) = Qa(0O, 001, 0O2)+Qa(0O2, 001, 0l)+<9a(01, 012, 0O2) + Qa(0O2, 012, 02) 

This can be expressed in terms of a matrix form as 

^?o(0O, 01, 02, 001, 002, 012) = T -Ra0 

The result of integrating out the new variables 0oi,0o2,0i2 is best under- 
stood by splitting the matrix R into a block L that acts on the original variables 
(0o,0i,02) another K acting on the new variables (0oi, 002 , 012) an d a third 
matrix J that mixes the two. 



R„ 



L J 
J T K 




Aa a -2a a 
K = I -2a a 4a a + 4/3 a 






-2a a 




"a 


Pa 





"a 




a a 





Pa 


a a 




The exponent in the above Gaussian can also be written as 

^o(0O, 01, 02, 001, 002, 012) = <tf> L(j) e + 20f J0i + (j>J Kfc 

(pi and e denote the set of internal (new) and external (original) vertices of 
the triangle. On completion of squares we get 

R a = <g £0 e + (0f + 0jJ T J FT- 1 )^(0 i + K- 1 J<j> e ) - <f e J T K- x J<i> e 
The result of the integral over the internal vertices then turns out to be 

Q2a = L a - J^K~ x J a 

§(5a Q + 6/3 Q ) \{-a a -P a ) \{-a a -2p a ) 

\{-a a -fla) U a + /3 a \{-Ct a -P a ) 

U-a a -2p a ) \{-a a -p a ) |(5a a + 6/3 a ) 



G 



3.1 Fixed Points of the transformation 

The above result should be equal to the quadratic form 

a 2 a (00 - 4>lf + (01 - 4>2f + /?2a(0O - §lf = 4> T Q 2a (j> 

On comparing its form we get 

a 2a = 7,{a a +P a ), l32a = ^(a a + 2/3 a ) 

If we define the ratio x a — ^- we can write these relations as 

1 + 2x 



X 2a = f(x a ), f(x) = 



4(1 + a;) 



a 2a = X(x)a a , X(x) = -(l + x) 

We now seek a fixed point of the transformation x a i-> x 2a which is a 
quadratic equation with two roots 



f(x) = x 

X = ; 



The solution with negative x should be discarded as it does not give an 
integrable Gaussian. The positive root corresponds to the scaling relation 

a 2a = Aa a , A = - ^3 + Vt^j 

Thus we set 

v/5- 1 



X= 4 

If we make the ansatz a a = aT^a.^ we get 



or 



/i = ° v - ' = — ; " ~ °- 611516 

log 2 log 2 

In summary, we have a scale covariant solution to the reproducibility condi- 
tion 



7 



^a(0o,0i,02) = C exp 



a-"ai | (0o - 0i) 2 + (02 - 0i) 2 + ^V^(0o - 4>2? 



4 

We can also determine the normalization constant using the relation 



^32(^ + 2^) ° 2a 
but that is less interesting. 

3.2 Correlation function 

Choose two points x\, x 2 on the plane separated by some distance L = \xi — x 2 \. 
Draw an iso- right triangle with the line connecting these points as a hypotnuse. 
We then bisect the right angle to get two iso-right triangles, and subdivide them 
and so on. After a large number of subdivisions we will get triangles of very 
small sides. We then integrate out all vertices expect the original triangle to 
calculate the correlation function: 

r( , _ JUi #« 11/ Fa W/)) - 0(^)] 2 

The integrals can be evaluated using the reproducibility property above. In 
the end we will end up with an integral over the two vertices of the original 
large triangle whose side is L : 



J dfodfadfa [0i ~ <h? exp \-^a, { ((0 O - 4>i? + (02 - 0i) 2 ) + ^(0o - 2 ) 2 )}1 
G(x lt ,x 2 ) = — ! 



fdfadfaexp [-^a, {((0 O - <M 2 + (0 2 - 0i) 2 ) + ^(0o " 2 ) 2 )} 
Thus, for some positive constant G, 



G(xi,x 2 ) = G\x\ - x 2 \^ ~ \xi - x 2 



10.611516 



Note that the magnitude of the correlations increase with distance. This is 
similar to the behavior of the average distance between two points on a Brow- 
nian motion. On the other hand this is not the behavior expected from a two 
dimensional quantum field theory with action S = — \ J \W<p\ 2 d 2 x = |(0, A(f>). 
The correlations behave logarithmically then 



/ 



(27T) 2 |fc| 2 81 1 21 



An interpretation is that it is a quantum field theory in a space of fractal 
dimension in between one and two. 
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Figure 4.1: 



4 Subdividing Tetrahedrons 

It is not possible to subdivide a tetrahedron into two that are similar to itself by 
a simple bisection: that would be equivalent to the ancient problem of doubling 
the cube which is impossible by ruler-compass methods (i.e., the cube root of 
2 cannot be found by such "Platonic" constructions.) But it is possible to 
subdivide a tetrahedron into eight tetrahedrons, each similar to the original. 
The trick is to choose right isosceles tetrahedra, whose faces are all right angled 
triangles, based on three mutually orthogonal edges of equal length|8 . The 



figure 4.1 shows an example of such a tetrahedron, the three mutually orthogonal 
edges being shown in bold lines. 
To be more specific, the vertices 

xq = (0,0,0), xi = (a, 0, 0), X2 = (a, a, 0), x 3 = (a,a,a) 



gives a right isosceles tetrahedron of smallest side a (Figure 4.2 ). Every face 



is a right angled triangle. There are three kinds of edges, of squared lengths 
a 2 , 2a 2 , 3a 2 . In fact, the squared length of the edge is just \i — j\a 2 , for i,j = 
0,1,2,3. To each such tetrahedron we associate a function T a (0O) (pi, 02, 03) 
which represents the result of integrating out the field in its interior. 

The internal vertices of the subdivision are the six points Xij = \ (xi + Xj ) 
for < i < j < 3. 

The reproducibility condition is the integral equation 



T2a{<P(i,(pl,(p2,<p3,) = J T a ((p ,(p 01 ,(p 2,(p03)T a ((pi,(p 01 ,(p 2,(p03) 

T a ((pl,<Pl2, 002, (p03)T a ((pl,(pl2, 013, <p03)T a ((p2, 012 , 002 , 003) 
T a ((p2, 012 , 013 , 003)^(02, 023, 013, (paz)T a {<Pz, 023, 013, 003) 
^001^002^003 ^012^013^023 
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Figure 4.2: 



Remark 1. It is important to have the correct order of arguments in the T a as 
we are not dealing with a regular tetrahedron. We have verified that with the 
orders given above, (different in detail from the reference above) the distance 
matrices are exactly half the distance matrix of the original tetrahedron. 

We assume a translation symmetry in field space (a global symmetry we 
expect for massless fields) 

T a {<t>a +4>,4>i + <f>> 02 + 0, 03 + 4>) = T a ((t>o, 4>i, 02, h) 
Also, assume a Gaussian ansatz for T a . 

T a (0 o , 0i, <j>2,fo) = e -Q«W°'^.« 

Qa(<h,<h.,<h,<h) = a((0 o -0i) 2 + (0 3 - 2 ) 2 ) +/3((0O-0 2 ) 2 + (03-0l) 2 ) 
+ 7 (0 O - 3 ) 2 + 5(0i - 2 ) 2 



This builds in the translation invariance as well as the reflection symmetry 
(0123) <=± (3210) of the right isosceles tetrahedron. 

The exponent of the integral (4.1) is a quadratic function of ten variables: 



Ra (00 1 01 1 02 , 03 1 001 > 002 , 003 , 012 , 013 , 023) = Qa (00 , 001 , 002 , 003 ) + Qa (01 , 001 , 002 , 003 ) 

+ Q a (01, 012, 002, 003) + Qa(01, 012, 013, 003) 
+ Qa(02, 012, 002, 003) + Qa(02, 012, 013, 003) 
+<5a(02,023,013,0O3) + Qa (03 , 023 , 013, 003) 
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This is equal to ^<f) T R(f> for some 10 x 10 symmetric matrix R, which can 
be calculated by evaluating the second derivatives (Hessian), a straightforward 
procedure in Mathematica. This matrix can be split into blocks 

R= ( J T k) 

where K is a matrix in the six dimensional subspace of the internal variables 
(0oi, • • • ! ^23), L is a matrix in the four dimensional subspace of the external 
variables (0o, ■ ■ ■ , <fo) and J is a 4 x 6 matrix that mixes them. The result of 
evaluating the integral over the internal variables by completing the squares is 
again a Gaussian with the Hessian matrix 

L — J T K~ x J 

4.1 Fixed Points of the transformation 

Thus we can write the renormalization group transformation as 

Qia = L a - Ja K a 

On comparing with the form of Q a , the parameters a 2a • • • &2a will be deter- 
mined as some homogeneous rational functions of a a ■ ■ ■ <5 a .The ratios 

/?2q 72a _^2a_ 
Oi2a ' a 2a ' a 2a 

depend only on the ratios x — ^f-,y — — ^f-; i.e. a a and a 2a can 

be scaled out. Thus we get a map of the projective space P 3 to itself. Each 
fixed point (x*,y*, z*) of this map is a Gaussian that is mapped into itself under 
our transformation, except for an overall multiplication of the quadratic form 
by ^ (the "eigenvalue" A) evaluated at (x*,y*,z*). A sensible solution must 
have a positive eigenvalue and also the parameters (x* ,y* , z*) must define a 
positive quadratic form Q a . The equations for the fixed points can be reduced 
to a polynomial for x 

-2+22x+132a; 2 +380a; 3 +699a; 4 +188a; 5 -1639x 6 -2488x 7 -948a; 8 +594x 9 +608a; 10 +152a; 1 

The fixed point values for the other variables arc determined in terms of this 
root. 



11 



—————(-88214437948 - 533552485049a; - 1626406118465a; 2 - 3995440895007a; 3 
107606064063 

-4531514030431a; 4 + 6599814153205a; 5 + 18490138345312a; 6 + 10764883382910a; 7 
-3123153862430a; 8 - 5175977687624a; 9 - 1495324502248a; 10 ) 

—-—^-——(-16477061906 - 533552485049a; - 1626406118465a; 2 - 3995440895007a; 3 
35868688021 

-4531514030431a; 4 + 6599814153205a; 5 + 18490138345312a; 6 + 10764883382910a; 7 
-3123153862430a; 8 - 5175977687624a; 9 - 1495324502248a; 10 ) 



Of the eleven solutions, only five are real and among them exactly one leads 
to a positive quadratic form Q a 



x « 0.699787174329122, y w 0.182495561519144, z w 2.547486684557433 

The eigenvalue i.e. the ratio for this solution is = A « 0.7958942. 

As for the right angled triangles, the correlation function can now be cal- 
culated. We can choose any two vertices on the original tetrahedron for this 
purpose and find the correlation between them. The choice of vertices is not so 
important as the distance between them is very large compared to the sides of 
the tetrahedrons obtained by subdivision. The field variables in those vertices 
obtained by subdivision is integrated out. The correlations will scale as A^ 1 
when the distance is doubled. That is, 

G(a; 1 ,a; 2 )^|x 1 -a; 2 |^~|a; 1 -a; 2 | - 329351 



5 Additional Examples 

5.1 Gaussian Measures 

In both cases above we get a correlation function that increases with distance. 
Are these exponents universal? Is it possible to find another subdivision pro- 
cedure that will produce a correlation that decreases with distance? To answer 
these questions we worked out several additional examples of subdivision. It 
turns out that each subdivision of the plane (into squares, triangles of different 
shapes etc.) gives a different exponent for the correlation function. Also, in 
most cases (except two) we get correlation functions that grow with distance. 



5.1.1 Subdividing a Square 



We can also subdivide the plane into squares (Figure 5.1 1. To each square of 
side a we can associate a Gaussian weight with quadratic form 
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Figure 5.1: * 



Qa{(f>l,hA3,M = Oi a [(01-0 2 ) 2 + (02-03) 2 + (03-04) + (^4-</'l) 2 ] 
+fla [(01 - 3 ) 2 + (02 ' <Pif] 



We subdivide each square by connecting the centroid X5 to the midpoints 
of the sides, x i2 , x 23 ,x 34:7 x±\. If we then integrate out the five new vertices, we 
will get the transformations as 

a 2a = -(2a a + a ), p 2a = Saa + A p a 

Fixed Points and Correlation Function We get the condition for a 
fixed point of by setting f(x a ) = x 2a where x a = 

A(2 + Ax + x 2 ) _ 
(2 + x)(8 + 4x) =X 

Of the three solutions, only x = — 1, \/3 — 1 lead to positive eigenvalues. 

a 2a = \a ai A = i,or ^ (l + V3j 
These two fixed points both give growing correlations 

G(xi - x 2 ) ~ \x-i - x 2 \^, n = 2, or 0.550016 
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5.1.2 Apollonian Subdivision 



Another example of a decreasing correlation function is the iterated Apollonian 
subdivision of a triangle: we join the vertices to an interior point, thereby 



subdividing the triangle into three triangles (Figure 5.2 1. Another interpretation 
is in terms of circles: we can regard the plane as being subdivided into three 
mutually tangent circles centered at the vertices of a triangle. Then we can 
insert a new circle in the interstitial region that is tangent to the other three 
circle. This process can then be repeated. Every tangency of a pair of circles 
corresponds to an edge connecting their centers. 

The scaling law of the correlation function can be computed just as before 
by computing the Gaussian integral over the new vertex x$ 



/ e -aa[(0o-0l) 2 + (0O-02) 2 + (01-02) 2 +2(0o-03) 2 +2(0 1 --03) 2 +2(02-03) 2 ]^ 3 

Each side emanating from vertex x% is shared by two triangles, which ex- 
plains the factor of two in those terms. The proportionality is given as the 
scaling is not done rigorously here. We say that one third of the vertices are 
integrated out roughly after each step and we get y/3 as the scale factor in two 
dimensions. 

Fixed Points and Correlation Function This leads to the transforma- 
tion 

5 
3< 
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and a correlation function with exponent 

G{x 1 ,x 2 ) ~ n = ~ -0.929947 

log V3 

Finally we see an example of a correlation decreasing with distance. This 
may be interpreted as a fractal with dimension greater than two. 

5.1.3 Dividing a triangle into six similar triangles 

Another example of a decreasing correlation function is the division of an equi- 



lateral triangle divided up into six equilateral triangles (Figure 5.3). Topolog- 
ically this is the same as the barycentric subdivision, or the bisection of the 
angles of the triangles 0: we add a new point at the centroid (on the in-center) 
and connect it to the original vertices and midpoints of the sides. But as long 
as the six smaller triangles occupy the same plane as the original one, it is not 
possible for them all to be equilateral. But we can assign equal lengths to all 
the edges as long as we do not insist that the graph has an isometric embedding 
in Euclidean space. By iterating the subdivision process, we will get a fractal 
which makes sense as a metric space, but not as a sub-manifold of Euclidean 
space. 

To each of the equilateral triangle we associate the Gaussian weight 
Q a (<j) , fafa) = a a [(<fo - <M 2 + (0i - <fe) 2 + (<fe - M 2 ] 

The field variables at the vertices (xoi, X02, X12) and the center of the original 
triangle (2:3) are integrated over to give a bigger equilateral triangle with the 
vertices (xq, xi, x 2 ) ■ A side of the newly formed triangle is bigger than the 
original triangle's side by a factor of v6- 

e-QvBaWoA,^) = f e -^(*o,0i,02,0 3 ,0oi,0o2,0i2) d 3d oid i2d i2 
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with 



Ra(<f>0, 01,02, 03, 001, 002, 012) = <3a(0O, 001, 03 ) + Qa(0O, 03, 002) + Qa(01, 03, 001 ) 

+ Qa(01, 012,03) + Qa(02,03,012) + Qa (02 , 002 , 03 ) 



Fixed Points and Correlation Function On computing the Gaussian 
integral and comparing it with the Gaussian form expected we get the following 
transformation 



a V6a = ~ a ° 



t 

The new quadratic form for the triangle formed will be 

Qy6a(0O, 01, 02) = \&a [(01 ~ O )' + (00 ~ 2 )' + (01 - <hf] 

Thus the correlation function scales with exponent as 

lo 4 

G(x 1 - x 2 ) ~ \x x - x 2 \> i , j-i = ° gl w -0.249078 

log V6 

Thus it is possible to get a correlation function that decreases with distance 
as well. 

5.2 Ising Models 

5.2.1 The One dimensional Ising Model 

We carry out the idea of sub division of the real line once more. However we 
are now interested in the Ising model now and our field values (o^) at the end 
points of the of the interval can be ±1. Instead of integrals we end up doing a 
discrete sum. Let the function associated with an interval be 



F a (a ,(Ti) = e «.(*o."i) 



There will not be any higher order terms of Oi since of = 1 and OiOj — 
icijk&k- We divide the interval into two halves now and the new variable is 
summed over. 



e Qa(q , 0.g-Ql) e Qa(q , 01.q , l) _ ^ e a o T0l(ff0+Tl)+27, 
(701 O 



2 7a 



e 
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Fixed Points and Correlation Function On comparing with the func- 
tion associated with the whole interval 

we get the following pair of equations considering all the combinations formed 
by the pair (ctq, cti) 



2e 2ja = e l2a ~ a2a 



e 27 " (e 2aa +e~ 2aa ) = e l2a+a2a 



The transformation obtained from these equations is 

p 2a a I — 2a n 



2 

The trivial fixed point a a — is obtained by setting a 2a = a a- The normal- 
ization constant can be determined from the relation between j a and ^y 2a . 

5.2.2 Ising Model on the Fractal Subdivision of a Right Angled Tri- 
angle 



We carry out the idea of sub division for a triangle now (Figure 5.4). Let the 
function associated with the triangle be 



F a (a ,a 1 ,a 2 ) = e Q »^°^ 

Qa(CTo,cri,Cr 2 ) = a a (a (Ji + <Ji<7 2 ) + Pai^O^) + la 



We now divide the original triangle in two new similar triangles. The new 
vertex is summed over and the total contribution from the two triangles is 
computed. 



^""^ e Qq(g'l,g'02,q-o) e Q a (<72:0-02,0-l) _ ^ gO„ 



! CT02 (ff +CT2 +2(Ti )+P a (Jl ((To +<T 2 ) + 2~^ a 

(T02 &02 



|^ e Qa(o"o+o-2+2cri) _|_ e -a a {a +a 2 +2a 1 ) 



e /3 a o-i(cr +c72)+27 a 



Fixed Points and Correlation Function On comparing it with the 
form of the energy associated with the bigger triangle 

F^ a {<J ,a U a 2 ) = e a ^a^(^+<T2) + ^ Q ( CT 0<T2)+7y2a 
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we get the following pairs of equations 



e 2 7a ^2a 0+e -2a a ) 

2e 2(7a-^a) 



In order to obtain these equations we look into all the tuples formed by as- 
signing values to (erg, <J\, <7 2 ) ■ Out of the set of eight equations we get only three 
independent equations. The transformations obtained from these equations are 



2 

= 2 ( e4 " a+e ' 4 "°) 

( e 2a a +e -2a a )2 



The trivial fixed points a a — 0, /3 = are obtained by setting = 
a a , /3y2a = Pa- 7a and 7y5a determine the normalization factor. 

5.2.3 Ising Model on sixfold subdivision of the triangle 

After having found trivial fixed points a a = till now we take up the example 
of the sixfold subdivision of the triangle to see if it has a non trivial fixed 
point. And indeed it has one. To each of the equilateral triangles we associate 
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a function as shown below. 

A sum is now carried over the vertices (<7oi, 002, 1712 > ^° f° rm a bigger 
equilateral triangle with the vertices (ctoj fi, 02)- 



— gQo(o-0,O-l,O- 2 ) 

= a a (cr cri + ctiCT 2 + a-2fo) +7o 



^ ' ^ ' ~y ' ^ ' e Qa(CTl,CTl2,CT3)gQa(ff2,CTl2,CT3) e Qa(CT2,CT02,CT3) e Q a (o-0,O-02.CT3)gQa(CT0,CT01,(T3) e Q a (cri,Cr01,O"3) 
O"01 cr 02 CT12 0"3 

e -6a a +6 7a ^2 _|_ e -2(-2+tT 2 +cro)aa _|_ g2(2+cr 2 +o- )a a _|_ e -2{-2+a 2 +a 1 )a a + e 2(2+a 2 +a 1 )a a 
^g-2(-2+cr +(Ti)Q a _|_ g 2(2+(To+o-i)Qa _|_ g-4(-3+CT2+o"o+o"i)a a ^ g4(3+(T2+crn+o'i)Q a 
_|_g-2(-4+2tT 2 +<To+cri)aa i g2(4+2tT 2 +o-o+(Ti)a a + e -2{-4+a 2 +2a +a 1 )a a + e 2(4+a 2 +2a +cr 1 )a a 
_|_ e -2(-4+CT 2 +<r +2o-i)a o _|_ g2(4+<j 2 +(j +2o-i)a a \ 



Fixed Points and Correlation Function The function associated with 
the bigger triangle should be of the form 

We get the following pairs of equations on comparing the two results. 

e -«VBa+TVB. = g-6(«a-7o) (3 + 6 g4a a + ^a a + 2( ,12a a + glGa^ 
e 3«VBa+TVB. = g-6(a a - 7a ) + g e 8a a + 3^6^ + g24 Qa ) 

Thus the renormalization group transformation on a is the iteration of the 
function 

3 + 6e 4a » + 4e 8Q » + 2e 12a » + e 16a - ' 
9 + 3e 8 «» + 3e 16Q - + e 24a - _ 

There is a trivial fixed point for the equation a a = f{a a ) at a a = and 
a non trivial fixed point at a a * — 0.159579. The normalization constants are 
determined by j a and 7 % /g a - 

To find the critical exponents we linearize the recursion relation at the fixed 
point [5]. 

/'(a„) = 1.93411 



/(") = -^-log 



19 



This shows that the fixed point is unstable. The scaling factor for this trans- 
formation is v6- Using these numbers we can calculate the thermal exponent 
y t — '"^ogjygj 1 ^ — 0.736312. This leads to the critical exponents for the correla- 
tion length and specific heat to be v = ~ = 1.35812, a = 2(1 — v) = —0.716239. 
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